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Underdoped cuprates exhibit a normal-state pseudogap, and their spins and doped carriers tend 
to spatially separate into 1- or 2-D stripes. Some view these as central to superconductivity, others 
as peripheral and merely competing. Using L^-rrSr^Cui-yZnyCU we show that an oxygen isotope 
effect in T c and in the superfiuid density can be used to distinguish between the roles of stripes and 
pseudogap and also to detect the presence of impurity scattering. We conclude that stripes and 
pseudogap are distinct, and both compete and coexist with superconductivity. 
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High-T c superconductors (HTS) remain a puzzle. Var- 
ious correlated states have been identified in HTS in- 
cluding antiferromagnetism, the pseudogap 1], nanoscale 
spin-charge stripesQ and, of course, superconductivity 
(SC). (Here we generalise "stripes" to include possible 
2D checkerboard structures |3|). The pseudogap is a 
nodal energy gap of uncertain origin that appears in the 
normal-state (NS) density of states (DOS). Its effects can 
be observed in many physical properties[l|,|^| . Several op- 
posing views are still current. One is that stripes play 
a central roleja], forming the pseudogap correlation [|j 
and/or mediating the SC pairing. Another is that the NS 
pseudogap arises from incoherent superconducting fluc- 
tuations which set in well above T C [?J. Another is that 
these states are independently competing^. Here stripes 
and pseudogap play a secondary role and SC is mediated 
by some other pairing boson. An unambiguous test of 
these opposing views is urgently needed. We show here 
that isotope effects provide such a test. 

The isotope exponent a(E) in a given property E is de- 
fined as a(E) = - (AE/E) /(AM/M), where M is the 
isotopic mass and E may be T c , the SC gap parameter, 
Ao, the pseudogap energy scale, E g , or the superfiuid 
density p s = A~ b 2 = p e 2 (n s / m* b ) . (X ab is the in-plane 
London penetration depth, n s is the carrier density and 
m* b is the effective electronic mass for in-plane trans- 
port). An isotope effect on T c was first discovered in 1950 
by Allen et al. for Sn 9]. They found a{T c ) w 0.5 ± 0.05 
which provided the central clue for the role of phonons in 
pairing and led 7 years later to the BCS theory of SC|10|. 

The situation with HTS is more complex. The oxy- 
gen isotope effect on T c was 

found[IH 

to be small, with 

a(T c ) « 0.06. However, with decreasing doping the effect 
rises and eventually diverges as T c — > 0[lj,[R|. Surpris- 
ingly, an isotope effect was also found in the superfiuid 
density^] (and attempts were made to resolve this into a 
dominant isotope effect just in m*0,0|). We will show 
that both of these unusual effects can be understood in 
terms of a normal-state pseudogap which competes with 



sciuj. We also predict and confirm an isotope effect 
in p s induced by impurity scattering. The isotope ef- 
fects in T c and p s are mapped as a function of doping in 
La2- 2; Sr a ;Cui_ a Znj,04 and we observe a canonical pseu- 
dogap behavior as well as a huge anomalous effect asso- 
ciated with stripes. The clear distinction between these 
effects shows that the pseudogap and stripe states are 
distinct and both compete with SC. 

An isotope effect, a(p s ), in the superfiuid density is 
surprising because for a simple BCS superconductor it is 
rigorously zero. According to Leggetts theorem, p s is just 
the total integrated spectral weight of the free carriers i.e. 
the total carrier density divided by the effective mass[is|. 
But, when there are strong departures from nearly- free- 
electron theory this need not be so. We identify two cases 
for HTS in which an isotope effect in p s arises: in the 
presence of (i) impurity scattering, and (ii) a pseudogap. 

HTS possess a d-wave order parameter and in the pres- 
ence of impurity scattering both T c and p s are diminished. 
The degree to which they are reduced depends upon the 
magnitude of the scattering rate, T, relative to the maxi- 
mum gap parameter, Ao, near k = (jr, 0). In the presence 
of a competing pseudogap, spectral weight removed by 
the pseudogap is no longer available for the condensate 
and, again, both T c and p s are diminished. The degree to 
which they are reduced depends upon the relative mag- 
nitudes of the pseudogap and the SC gap. Thus, T c and 
p s are reduced according to the magnitude of the ratios 
r/A for impurity scattering, and E g /A for a pseudo- 
gap. In either case, a relatively small isotope effect in 
Ao will necessarily produce enhanced isotope effects in 
T c and p s which diverge as T c — > 0. Now it has been 
shown from specific heat, NMR and ARPES that, with 
increasing doping, E g decreases and closes abruptly at a 
critical doping state, p CT it — 0.19 holes/Cu, in the lightly 
overdoped regime^. It follows that the isotope effect in 
Ps should disappear at critical doping where the pseu- 
dogap closes provided that impurity scattering is absent. 
Our initial task is to quantify these effects. 



2 



C = r/r 

r. 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



0.0 0.2 0.4 0.6 



1.0 





1.0 


o 


8 


K° 


0.6 






o 




o 


0.4 


Cl 




CL 


0.2 









1.0 


o 

h-° 


0.8 




0.6 






o 






0.4 


o 


ex 






0.2 




0.0 







(a) 




-v TIT 

C CO 




P 'P 




- 


1 r- 


1 h- 

. TIT" 


(b) ■ 


p /p X 













0.0 0.2 0.4 0.6 0.8 

? = E/(2.397kX 



1.0 



FIG. 1: The fractional suppression of T c and superfluid den- 
sity in a d-wave SC with (a) unitary-limit impurity scattering 
and (b) a NS triangular pseudogap with gap energy E g . 



First, we recall that there is no isotope effect in the 
pseudogap. We have examined the 89 Y Knight shift in 
YBa2Cii408 using magic angle spinning with extremely 
narrow linewidths (sslOO Hz) and found no isotope effect 
within the bounds a(E g ) < 0.01 0|- Though not essen- 
tial, we proceed under the assumption that an isotope 
effect is confined to the pairing gap, Ao, and absent from 
the pseudogap, E g . The small isotope effect observed in 
X/T\T in the same compound does not reflect an isotope 
effect in the pseudogap. Using the enhanced suscepti- 
bility formalism it devolves, rather surprisingly, into an 
isotope effect in the paramagnon frequency 

Theory. Impurity scattering for a d-wave order param- 
eter has been investigated by many authors [13. l20| . We 
summarize the results in Fig. 1(a) which shows the de- 
pression of T c and p s as a function of £ = r/r c , where 
r c is the critical scattering rate for fully suppressing 
T c . The reduction in T c follows the standard Abrikosov- 
Gorkov equation. In the unitary limit, the scattering rate 
r = ni/irN(Ep) where rij is the density of scatterers and 
N(E) is the NS DOS. Fig. 1(a) shows that T C (C) /T c0 falls 
at first slowly then accelerates while p s (() falls at first 
rapidly then slows as T grows. We define the functions h 
and g given by p s (()/p s o = h(C) and T c (()/T c0 = g((). 
The isotope effects in p s and T c are 



a( Ps ) = -C(h'/h)a(T c0 ) 
a(T c ) = [l-((g'/g)]a(T c0 ) 
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FIG. 2: Oxygen isotope exponent in superfluid density plot- 
ted against that in T c for I^-nSr^Cui-yZnyO^. Panel (a) 
shows the calculated (red line) and observed effect of impuri- 
ties for overdoped (x = 0.19, y = 0, 1, 2, 2.4%) and strongly 
underdoped (x — 0.09, y = 0, 0.5, 1%) samples. Previous 
data is also shown (blue down triangles). Panel (b) shows 
more data for 0.10 < x < 0.19, with y — (blue squares), 
and 1% Zn (green squares) and reveals an anomalous devia- 
tion from the canonical pseudogap line (red line). Previous 
data is shown for Yi_ z Pr z Ba2Cus07-5 (blue up triangle). 



The prime indicates differentiation of or g(C)- Thus 



a(p s ) = -C(/ l 7^)[l-C(.97.9r 1 «(r c ) 



(2) 



We have shown previously 0] that the isotope ef- 
fect in T c across the entire phase diagram is consistent 
with an underlying exponent (in the absence of scatter- 
ing and pseudogap) of a(T c o) w 0.06. The red line in Fig. 
2(a) shows a(p s ) plotted versus a(T c ) using this value. In 
the absence of impurity scattering a(p s ) = and a(T c ) 
= a(T c o) — 0.06. This is the left-hand termination of 
the red line. With increasing scattering both a(p s ) and 
a(T c ) rise along the line, and finally diverge as T — > T c . 

Turning to the pseudogap, specific heat|8( and tunnel- 
ing measurements |2l] show that the pseudogap is non- 
states-conserving, with an approximately triangular en- 
ergy dependence, and pinned to the Fermi level, Ep. We 
assume therefore a triangular normal-state DOS: 



N(E) 



N x 
N ; 



\E-E F \/E g (p); 



\E- 
\E 



E F \<E g (p), 
E F \>E g (p). 



(3) 

and solve standard weak-coupling d-wave BCS expres- 
(1) sions to calculate T c as a function of E g . For this par- 
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ticular NS DOS T c -> as E g -> 2.397k B T° where k B 
is Boltzmanns constant and T c ° = T c (E g — 0). Fig. 1(b) 
shows T c plotted as a function of ( = E g / (2.397k B T°). 
As for impurity scattering, the depression in T c is slow 
at first and more rapid as £ — ► 1. 

Elsewhere [j^ we have calculated the effect of a trian- 
gular pseudogap on p s . The approach is admittedly for 
a Fermi liquid but we note that the effects we describe 
are dominated by the nodal regions of the Fermi surface 
where such a Fermi liquid approach is more likely to be 
valid. p s (£) is plotted as a function of £ in Fig. 1(b). This 
exhibits an initial rapid fall which slows as E g grows and 
Q — > 1. The overall behavior is qualitatively similar to 
that shown in Fig. 1(a) for impurity scattering, but dif- 
fers in detail. We could therefore define new functions 
/i(£) and <?(£) as above and derive an equation formally 
identical to eq. (2) to describe the isotope effects in p s 
and T c associated with the presence of the pseudogap. 
These equations show that when the pseudogap closes 
at critical doping we have £ = and a(p s ) — while 
a(T c ) = a(T®). The resultant curve a(p s ) versus a(T c ) 
almost exactly coincides with the red line in Fig. 2(a). If 
there were an isotope effect in the pseudogap a(T c o) in 
eq. (1) should be replaced by [a(T c0 ) — a(E g )} and eq. 
(2) and the red line in Fig. 2 remain unchanged. 

Experimental details. La2- a; Sr :l; Cui_yZn, / 04 samples 
were synthesized by solid state reaction at 985°C in air by 
repeated milling, pelletization and reaction until phase 
pure as determined by x-ray diffraction. Two small, ap- 
proximately 2x2x3 mm 3 bars, were cut from alongside 
each other at the centre of each of the resultant pellets to 
ensure, as much as possible, identical pairs. They were 
isotope exchanged in identical quartz tubes, one charged 
with 16 and the other with 18 0, side by side in a furnace. 
The 18 O gas (from Isotec) was 99% enriched and several 
exchanges were employed until about 95% exchange was 
achieved. On the final exchange the samples were slow 
cooled then annealed for 15 hours at 500°C to ensure oxy- 
genation to full stoichiometry. The degree of exchange 
was confirmed by Raman measurements of the spectral 
shifts of the oxygen phonons. 

To determine the isotope shifts in T c and p s we car- 
ried out field-cooled DC magnetization measurements in 
the mixed state at 150 Oe. For this regime Zhao and 
Morris ^ adopted the relation |23j| 

(-M) 1 / 2 oc (r g /X)[T c (H)-T]x[\dH c2 /dT\/(2 K 2 ~l)T c ] 1/ ' 

(4) 

for the limit, near T Cl of A 3> r g . Here r g is the mean 
radius of the SC grains. These authors showed that this 
relation could be used to deduce separate isotope effects 
in n s and to*. But the algebra was incorrect (see Ap- 
pendix below). A further problem arises [2^ in that, for 
small particles, this relation does not satisfy the sum rule, 
/io / M(H)dH = Uo = the condensation energy. With a 
mean grain size of 25/im and A a b(0) ranging from 0.2 to 
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FIG. 3: Magnetization versus T for La2-xSr :r Cui_j / Znj / 04. 
Panel (a) shows x = 0.19 with y = 0, 1, 2, 2.4% (solid curves) 
and x = 0.09 with y = 0, 0.5, 1% (dashed curves). Panel (b) 
shows data for Zn-free samples with 0.09 < x < 0.22. 

0.32/im 25], we adopt the limit A <C r g which is clearly 
satisfied up to a few K below T c . This yields 

-Mac \~ 2 [l-T/T c \. (5) 

Thus the isotope coefficient in the slope of —M(T) is 
given by a(p s ) — a(T c ). What we report is the partial 
isotope exponent due to the change in oxygen mass only. 
We obtained qualitatively similar results with Meissner 
state measurements at 10 Oe (not shown). 

Results. We start first with the effect of impurity scat- 
tering in the overdoped region x — 0.19 where the pseu- 
dogap is absent. Illustrative plots of magnetization ver- 
sus temperature are shown in Fig. 3(a) for y = 0, 1, 2 
and 2.4% (solid curves). It evident that an isotope ef- 
fect in T c is present in each but that one in p s is only 
present in the higher Zn concentrations. Values of a(p s ) 
are plotted versus a(T c ) in Fig. 2 (a) (left-hand cluster 
of blue, green, black and mauve squares) and they are 
seen to be roughly consistent with the model calculation. 
The fact that a(p s ) — > as y — ► indirectly shows that 
any disorder potential present in Zn-free La2- a; Sr 2 ;Cu04 
in the overdoped region is too small to present significant 
scattering and hence a(p s ) = 0. It also seems unlikely 
that there is any significant phase separation because the 
domain walls would surely act as scattering centers. 

Turning to the heavily underdoped pseudogap region, 
we show magnetization curves in Fig. 3(a) for x = 0.09 
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with y = 0, 0.5 and 1% Zn (dashed curves). The resul- 
tant a(p s ) versus a(T c ) values are plotted in Fig. 2(a), 
shown by the blue, red and green data points to the 
right. These continue to track up the canonical curve, 
showing that the pseudogap and impurity scattering have 
essentially the same effect in such a plot. To this data 
we add previously-reported [26j values for I^-zSr^CuC^ 
obtained using muon spin relaxation (/iSR) with x — 
0.080 and 0.086 (blue down triangles). The collective 
data is generally consistent with the model. 

Fig. 3(b) shows a selection of illustrative plots of M vs 
T for Zn-free samples with x ranging from 0.09 to 0.22. It 
is immediately evident from the low-T values of M that 
a(p s ) = for all x > 0.19 but becomes non-zero and 
large as x falls below 0.19. Values of a(p s ) are plotted 
against a(T c ) in Fig. 2 (b) (blue squares) and increas- 
ing doping is shown by the arrow. Here, a remarkable 
anomaly is evident. The overdoped data and the heavily 
underdoped data lie near the canonical pseudogap line. 
But near x = 0.12 the data deviates drastically from this 
canonical behavior. This is presumably due to the pres- 
ence of charged stripes, inferred from neutron scattering 
near p=l/8, which provide strong electronic coupling to 
the lattice. If the pseudogap itself arose from fluctuat- 
ing stripes one might expect the anomaly to drive up the 
canonical line. The huge deviation suggests a fundamen- 
tally different behavior and clearly distinguishes stripes 
from the pseudogap near p=l/8. 

In order to further test this interpretation we examined 
the effects of non-magnetic Zn substitution. Our expec- 
tation was that the combined effects of spin vacancies 
and the tendency of Zn to enhance the canonical behav- 
ior would be to broaden and weaken the anomaly pushing 
it up the canonical line. Fig. 2(b) shows the effect of 1% 
Zn substitution (green squares). The contour, indicated 
by the green curve, confirms our expectations. 

Finally, we show by the upward open triangle in 
Fig. 2(b) recently reported a(p s ) and a(T c ) data for 
Yi_zPr z Ba2Cu307_5 obtained by Khasanov et al. using 
uSR|27j. For z = 0.3 these authors found T c = 59.3K 
and we estimate that the doping state is very close to 
p = 0.125. And yet the data resides close to the canon- 
ical line completely free of the anomalous deviation as- 
sociated with stripes. It is clear from inelastic neutron 
scattering studies that the YBa2Cu 3 07_5 compound ex- 
hibits a much weaker tendency to stripe formation. Con- 
sistent with this we find this sample exhibits essentially 
stripe-free canonical pseudogap behavior. 

We conclude that our results and analysis demonstrate 
a clear distinction between the canonical effects on the 
superfluid density arising from the pseudogap and impu- 
rity scattering on the one hand and stripe correlations 
on the other. We achieve this by examining a plot of 
a(p s ) versus a(T c ) which is relatively insensitive to the 
precise details of the NS DOS. Stripes cause a huge de- 
viation from this canonical behavior associated with the 



strong electronic coupling to the lattice arising from spa- 
tial charge modulation. On the basis of these results 
we make the important conclusion that stripe and pseu- 
dogap correlations are fundamentally different and both 
compete with each other and with superconductivity. 

We acknowledge financial support from the Mars- 
den Fund and the MacDiarmid Institute (JLT, JS and 
GVMW) and from Trinity College, Cambridge and the 
Cambridge Commonwealth Trust (RSI). 

Appendix - isotope effect in m*? 

Several authors have considered the possibility 

that the isotope effect in p s oc n s /m* may be resolved 
into a(n s ) - a(m*) and they have sought to determine 
these two components separately. 

Zhao et a/.[lfj investigated the oxygen isotope depen- 
dence of the orthorhombic/tetragonal (O/T) transition 
in La2- a; Sr a ;Cu04 and found a null effect. Because the 
O/T transition temperature is doping dependent they 
took this to indicate that there was no isotope effect in 
the carrier concentration and consequently the isotope 
effect in p s derives wholly from the isotope effect in to* 
i.e. at(ps) = —a(m*). However, the location of the O/T 
transition is an ion-size dependent effect not primarily a 
doping effect and, moreover, there is no simple relation- 
ship between the doped hole concentration, x, and the 
carrier concentration. 

Elsewhere, Zhao and Morris^jJ use eq. (5) for 10G 
measurements to yield a magnetisation slope 

Pi cx r 2 g n s /T c m* , (6) 

while for 150G measurements they use eq. (4) from which 
they deduce 

P 2 <xr g n 5 J 3 /T c m*. (7) 

Clearly, measurement of the isotope effects in Pi and P 2 
would allow extraction of the individual isotope effects in 
n s and in to*. However, eq. (4) does not lead to eq. (7). 
To see this we consider the relation [isj 

H c2 (0) = 0.7 T c [dH c2 /dT] Tc = O /(2<(O) 2 ). (8) 

which, on substitution in cq.(4) when k = A(0)/£(0) ^> 1, 
reduces to eq. (5) and 

Pi oc Pi oc n s /T c m* . (9) 

Thus the isotope effect in p s cannot be separated into 
separate contributions from a(n s ) and a (to*) in the way 
suggested by Zhao and Morris. 
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